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Oblique  Propagation  of  Groundwaves  Across  a  Coastline — 

Part  III 


James  R.  Wait 


Contribution  from  the  Central  Radio  Propagation  Laboratory,  National  Bureau  of  Standards,  Boulder,  Colo. 

(Received  November  8,  1963) 

This  paper,  w  hich  is  a  continuation  of  two  earlier  papers  of  the  same  title,  contains 
numerical  results  for  the  field  anomaly  hear  a  coastline  when  the  surface  impedance  changes 
•in  a  linear  manner  between  land  and  sea.  The  earlier  results  for  an  abrupt  boundary  are 
recovered  as  the  width  of  the  transition  region  is  reduced  to  zero.  Iii  general,  it  is  found 
that  the  characteristics  of  the  transition  region  will  not  produce  significant  modifications  of 
the  transmitted  field.  However,  the  magnitude  of  the  reflected  field  is  greatly  reduced  as 
the  Width  of  the  transition  zone  is  increased  beyond  about  one-quarter  wavelength. 


1.  Introduction 

lit  part  1  [Wait,  1 963]  of  a  series  of  papers  of  this 
title,  the  propagation  of  radio  waves  across  a  flat- 
lying  coastline  was  investigated  theoretically.  In 
part  II  [Wait  and  Jackson,  1963];,  the  influence  of  an 
elevation  change  between  land  and  sea  was  con¬ 
sidered.  It  is  the  purpose  of  the  present  paper  to 
consider  again  the  flat-lying  coastline  but  with  special 
attention  given  to  a  nonabrupt  variation  of  the 
conductivity  at  the  coastline. 

As  indicated  in  part  1,  the  assumption  of  a  sudden 
change  of  electrical  properties  gives  rise  to  an  appar¬ 
ent  singularity  of  the  field  right  at  the  coastline.  It 
is  shown  in  what  follows  that  the  “singularity”  is  not 
present  when  the  surface  impedance  changes  grad¬ 
ually  at  the  junction  or  coastline. 

2.  Formulation 

The  model  is  the  same  as  that  used  in  part  I.  The 
essential  features  are  described  briefly  here.  The 
transmitting  antenna,  A,  located  on  the  land,  is  at  a 
great  distance  from  the  coastline.  For  present  pur¬ 
poses,  it  is  only  necessary  to  assume  that  A  is  a  source 
of  vertically  polarized  groundwaves.  The  receiving 
antenna,  B,  is  relatively  near  the  coastline  hut  it  may 
be  on  the  sea.  It  is  further  assumed  that  B  is  equiva¬ 
lent  to  an  infinitesimal  vertical  dipole  and  thus  it 
responds  only  to  the  vertical  electric  field. 

From  consideration  of  reciprocity,  it  is  dear  that 
the  role  of  transmitter  awl  receiver  may  be  inter¬ 
changed.  Thus,  in  general,  it  is  most  meaningful  to 
express  the  results  in  terms  of  mutual  impedance 
zm+Azm  between  the  respective  terminal  pairs  of 
antennas  A  and  B.  Here  zm  is  the  mutual  impedance 
if  the  surface  impedance  were  a  constant  Z  for  all 
points  of  the  earth’s  surface.  Thus,  Az„,  is  the 
modification  of  the  mutual  impedance  which  results 
from  the  presence  of  the  inhoiuogeneily  of  surface 
impedance  Z'.  In  the  present  problem,  the  surface 


impedance  of  the  land  is  Z  while  that  of  the  sea  is 
Z' .  The  situation  is  illustrated  in  figure  1  where  the 
earth’s  surface  is  the  (x  y)  plane  of  a  Cartesian  coor¬ 
dinate  system  and  the  coastline  is  at  x=0. 

3.  Statement  of  Formulas 

Quoting  from  part  I,  the  working  formula  for  the 
mutual  impedance  increment  Azm  is  given  by 


XimkCWx-fiWdx,  (1) 

where 


fix)  =  (Z'  —Z)fi„,  7/0  =  120ir, 

II  '[f  is  the  Hankcl  function  of  order  zero  of  the 
second  kind, 

C,=cos  0,„  A  =2*7  wavelength, 

shortest  distance  from  coastline  to  B  (positive  if 
B  is  over  sea  and  negative  if  B  is  over  land) , 

and,  finally,  e  is  chosen  to  be  sufficiently  large  that 
fix)  art)  for  xfi  —  e. 

For  purposes  of  the  present  paper,  some  further 
simplifications  may  lie  made.  Specifically,  it  is 
noted  that 

f(x)  =  —A<i(x)ci,l\  (2) 

where  A0(x)  is  approximately  a  real  positive  quantity 
whose  magnitude  is  small  compared  with  unity  (e.g., 
A(i<^0.1).  For  example,  at  sufficiently  large  positive 
values  of  x, 
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finite  transition  region  of  width  </0  where  Ao(z)  Varies 
in  a  linear  manner.  Thus 

Ao(z)=0  for  £<— d(l/2, 

=  [?±^o/?)-]  Ao  for  — (/0/2<z<d0/2, 

— Ao  for  z>f/0/2.  (5) 

With  this  model,  the  surface  impedance  is  imagined 
to  change  in  a  linear  manner  over  the  distance  </0  as 
indicated  in  figure  lb. 

To  simplify  the  computational  problem,  certain 
dimensionless  quantities  are  introduced  as  follows. 

a—kC\X,  on=kCA,  D=kCid0, 

and  D0=DtC,—k<l0.  Then,  it  easily  follows  that  (1) 
may  be  written  in  the  form 


A(X) 


Az„  AqCQt/i  . 


-  “1 1 1'r " 


j  Htr'Wa^c 

(6) 


where,  in  the  case  of  a  linear  variation, 
G= 0  for  a<— 7)/2, 

=  for  -m<«<bl2, 


G(a) 


Fltil  RB  11),  e.  The  linear  .surface  impedance  variation. 


where  <r,  and  e„  are  the  conductivity  and  permit  tivity 
of  the  land  and  <r'  and  t'  are  the  conductivity  and 
permittivity  of  the  sea.  In  most  cases  of  practical 
interest,  u  is  sufficiently  low  and  a'  is  sufficiently 
large  that 

oj/«tk)i/2,  (4) 

which  is  real. 

4-  Linear  Variation  of  Surface  Impedance 

Ill  part  I,  it  was  assumed  that  A,iO)  was  a  step 
function  at  z=0  corresponding  to  an  abrupt  bound¬ 
ary.  Here  it  is  desirable  to  allow  A„(z)  to  have  a 


=  1  for  a>/)/2.  (7) 

Therefore, 

f  -;)2<„<)),2, 

=  0  for  «<  — 1)12  and  «^>/12,  IS) 

as  indicated  in  figure  le  where  G  is  plotted  versus  <*. 
I  The  mutual  impedance  formula  now  becomes 

A,,  [r^  +^  Aotf ,  (9) 

where 

Q,  2]J  (10) 

and 
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To  effect  these  integrations,  it  is  desirable  to  make 
use  of  the  following  two  identities: 

[xe±‘W (x)  ? \HP (x)l]=«±W  (x) ,  (12) 

and 

“  \xH^  (x) +(1  ?ix)ff,®  (x)]l=xe±frHo2’  (x) , 

(13) 

where  upper  (or  lower)  signs  are  to  be  considered 
together.  In  verifying  these  identities,  it  is  neces¬ 
sary  to  employ  two  well-known  relations  in  Bessel 
function  theory  [MacLachlan,  1934].  These  are 

=  (14) 

and 

^fl,{B(x)=fly>(x)“^^-  (15) 

It  is  now  apparent  that  A zm/zm  may  be  expressed 
entirely  in  terms  of  ttankel  functions  of  order  zero 
and  one,  with  various  real  arguments.  Using  stfbh 
expressions,  curves  of  the  real  and  imaginary  parts 
of  Q  have  heen  prepared.  Some  of  these  are  shown 
in  figures  2a  to  8b  where,  in  each  case,  the  abscissa 
is  the  parameter  a,  or  kC\(lx. 

For  purposes  of  discussion,  it  is  convenient  to  de¬ 
scribe  Q  as  the  field  anomaly  resulting  from  the 
inhomogeneity.  In  view  of  the  relation 

— ^=Ao  [Re  Q+i  Im  Q],  (16) 

Zm 

it  is  evident,  for  Ao  real,  that  Re  Q  is  a  measure  of  the 
amplitude  change  of  the  field  where  Iin  Q  is  a  meas¬ 
ure  of  the  phase.  While  the  results  are  still  valid  if 
Ao  is  complex,  the  above  simple  interpretation  is  not 
applicable.  Also,  it  should  be  remembered  that  under 
all  conditions,  \AzJ  2„|<<1  or  |Ao<?l<<l>if  the  pres¬ 
en  t  results  are  to  be  given  any  confidence. 

5.  Discussion  of  Numerical  Results 

The  curves  in  figures  2a  and  2b  are  applicable  to 
normal  incidence.  The  sinusoidal-like  ripples  for 
negative  values  of  at  may  be  regarded  as  an  interfer¬ 
ence  pattern  resulting  from  the  combination  of  the 
incident  wave  and  the  reflected  wave.  It  is  apparent 
that,  as  the  electrical  width  D„  of  the  transition  zone 
is  increased,  the  magnitude  of  the  reflected  wave  is 
generally  decreased.  On  the  other  hand,  the  char¬ 
acteristics  of  the  transmitted  wave,  for  large  positive 
values  of  a,,  are  not  appreciably  modified.  However, 
it  may  be  observed  that  in  the  proximity  of  the  coast¬ 
line,  the  nature  of  the  field  is  profoundly  influenced  by 
the  width  of  the  transition  region.  In  general,  the 
rapid  variations  of  the  field  are  smoothed  out  when 
the  transition  distance  Da  is  increased.  As  may  be 


Figure  2.  The  real  and  the  imaginary  parts  of  Q  as  a  function 
of  or,  or  KC|  d,  where  d,  is  the  distance  from  the  coastline. 


seen  in  figure  1 ,  the  field  is  beginning  to  have  a  singu¬ 
lar  behavior  for  the  smallest  value  of D0  (i.e.,  1)  which 
is  very  similar  to  the  corresponding  curves  given  in 
part  I  for  the  abrupt  boundary  (i.e.,  Z?0=  0). 

The  curves  in  figures  3a  and  3b  for  0„=2O°  are  very 
similar  to  those  in  figures  2a  and  2b  for  0o=O°.  Evi¬ 
dently,  a  slight  obliquity  in  the  traverse  across  the 
coastline  does  not  modify  noticeably  the  field  be¬ 
havior.  However,  for  more  oblique  conditions,  the 
curves  in  figures  4a,  4b,  5a,  and  5b  for  0o=4O°  and  50° 
indicate  that  significant  modifications  take  place  for 
angles  iti  this  range.  Principally,  it  is  noted  that 
when  0O  is  in  the  vicinity  of  459  the  reflected  wave  is 
greatly  reduced  in  magnitude.  This  effect  was  also 
present  in  the  case  of  an  abrupt  boundary. 

The  curves  in  figures  5a  to  Sb  illustrate  the  be¬ 
havior  of  the  field  anomaly  for  rather  highly  oblique 
angles.  Again,  as  in  the  case  of  the  abrupt  boundary, 
the  reflected  wave  is  quite  strong.  Furthermore,  the 
magnitude  of  the  reflected  wave  is  not  appreciably 
reduced  when  P„  is  increased  from  l  to  5.  However, 
in  nearly  all  angles  0„  the  behavior  of  the  transmitted 
field  is  not  appreciably  modified  by  changes  in  the 
width  of  the  transition  zone. 


Figure  3.  The  real  and  the  imaginary  ■parte  of  Q  os  a  function 
of  ai  or  KCi  di  where  d|  is  the  distance  from  the  coastline. 


The  relative  insensitivity  of  the  Q  versus  a,  curves, 
at  highly  oblique  angles,  to  variations  in  />„  is  reminis¬ 
cent  of  the  reflection  of  waves  from  horizontally 
stratified  media  [Wait,  1962],  In  the  latter  case  it  is 
known  that  the  influence  of  diffusiveness  or  non¬ 
sharpness  of  the  boundary  is  minimized  at  highly 
oblique  incidence. 

6.  Conclusions 

The  results  given  in  this  paper  would  seem  to  shed 
considerable  light  on  the  nature  of  the  vertical  elec¬ 
tric  field  near  the  boundary  of  separation.  The 
singularly  or  infinite  behavior  of  the  field  obtained 
in  earlier  studies,  is  not  present  when  the  surface 
impedance  changes  gradually  between  the  land  and 
sea  portions.  However,  as  the  transition  zone  is 
diminished  the  field  becomes  rapidly  varying  in  the 
vicinity  of  the  coastline  and  becomes  very  similar  to 
the  predicted  behavior  for  an  abrupt  boundary.  A 
conclusion  similar  to  this  was  arrived  at  bv  Godzjnski 
[19621  whose  analysis  communicated  to  me  in  outline 
was  restricted  to  normal  incidence  (i.e.,  0n=O°). 


294 


295 


